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Abstract. Let 5 C M, k+m be a compact semi-algebraic set defined by Pi > 

0, P t > 0, where P t G M[Xi, . . . ,X k ,Yi, . . . , Y m ], and deg(P,) < 2, 1 < i < 

1. Let it denote the standard projection from R* +m onto R m . We prove that 
for any q > 0, the sum of the first q Betti numbers of 7r(S) is bounded by 
(k + m)°( ql \ We also present an algorithm for computing the the first q Betti 
numbers of n(S), whose complexity is (k + to) 2 5 ' . For fixed q and £, both 
the bounds are polynomial in k + m. 



1. Introduction 

Designing efficient algorithms for computing the Betti numbers of semi-algebraic 
sets is one of the outstanding open problems in algorithmic semi-algebraic geometry. 
There has been some recent progress in this area. It has been known for a while 
that the zero-th Betti number (which is also the number of connected components) 
of semi- algebraic sets can be computed in single exponential time. Very recently, 
it has been shown that even the first Betti number, and more generally the first 
q Betti numbers for any fixed constant q, can be computed in single exponential 
time |l()l |H] . Since the problem of deciding whether a given semi-algebraic set in 
R fc is empty or not is NP-hard, and that of computing its zero-th Betti number 
is #P-hard, the existence of polynomial time algorithms for computing the Betti 
numbers is considered unlikely. 

One particularly interesting case is that of semi- algebraic sets defined by qua- 
dratic inequalities. The class of semi-algebraic sets defined by quadratic inequal- 
ities is the first interesting class of semi-algebraic sets after sets defined by linear 
inequalities, in which case the problem of computing topological information re- 
duces to linear programming for which (weakly) polynomial time algorithms are 
known. From the point of view of computational complexity, it is easy to see that 
the Boolean satisfiability problem can be posed as the problem of deciding whether 
a certain semi-algebraic set defined by quadratic inequalities is empty or not. Thus, 
deciding whether such a set is empty is clearly NP-hard and counting its number 
of connected components is #P-hard. However, semi-algebraic sets defined by qua- 
dratic inequalities are distinguished from arbitrary semi-algebraic sets in the sense 
that, if the number of inequalities is fixed, then the sum of their Betti numbers 
is bounded polynomially in the dimension. The following bound was proved by 
Barvinok pjj. 

Key words and phrases. Betti numbers, Quadratic inequalities, Semi-algebraic sets, Spectral 
sequences, Cohomological descent. 
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Theorem 1.1. Let S C E fe be a semi-algebraic set defined by the inequalities, 

k 

Pi > 0,-..,Pt > 0, deg(Pi) < 2,1 < i < I. Then, ^h(S) <k° {l \ where 

bi(S) denotes the i-th Betti number, which is the dimension of the i-th singular 
cohomology group of S, H Z (S;Q), with coefficients in Q. 

In view of Theorcm ll.il it is natural to consider the class of semi- algebraic sets 
defined by a fixed number of quadratic inequalities from a computational point of 
view. Algorithms for computing various topological properties of this class of semi- 
algebraic sets have been developed, starting from the work of Barvinok [2 , who 
described a polynomial time algorithm for testing emptiness of a set defined by a 
constant number of quadratic inequalities. This was later generalized and made 
constructive by Grigoriev and Pasechnik in where an algorithm is described 
for computing sample points in every connected component of a semi-algebraic set 
defined over a quadratic map. More recently polynomial time algorithms have 
been designed for computing the Euler-Poincare characteristic (Hj as well as all the 
Betti numbers |7] of sets defined by a fixed number of quadratic inequalities (with 
different dependence on the number of inequalities in the complexity bound) . Note 
also that the problem of deciding the emptiness of a set defined by a single quartic 
equation is already NP-hard and hence it is unlikely that there exists polynomial 
time algorithms for any of the above problems if the degree is allowed to be greater 
than two. 

A case of intermediate complexity between semi- algebraic sets defined by poly- 
nomials of higher degree and sets defined by a fixed number of quadratic sign 
conditions is obtained by considering projections of such sets. The operation of 
linear projection of semi-algebraic sets plays a very significant role in algorithmic 
semi-algebraic geometry. It is a consequence of the Tarski-Seidenberg principle 
(see for example [TT], page 61) that the image of a semi-algebraic set under a 
linear projection is semi-algebraic, and designing efficient algorithms for comput- 
ing properties of projections of semi-algebraic sets (such as its description by a 
quantifier-free formula) is a central problem of the area and is a very well-studied 
topic (see for example jHJ or 1 1 1| - Chapter 14). However, the complexities of the 
best algorithms for computing descriptions of projections of general semi- algebraic 
sets is singly exponential in the dimension and do not significantly improve when 
restricted to the class of semi-algebraic sets defined by a constant number of qua- 
dratic inequalities. Indeed, any semi-algebraic set can be realized as the projection 
of a set defined by quadratic inequalities, and it is not known whether quantifier 
elimination can be performed efficiently when the number of quadratic inequalities 
is kept constant. However, we show in this paper that, with a fixed number of 
inequalities, the projections of such sets are topologically simpler than projections 
of general semi-algebraic sets. This suggests, from the point of view of designing 
efficient (polynomial time) algorithms in semi-algebraic geometry, that projections 
of semi-algebraic sets defined by a constant number of quadratic inequalities is the 
next natural class of sets to consider, after sets defined by linear and (constant 
number of) quadratic inequalities, and this is what we proceed to do in this paper. 

In this paper, we describe a polynomial time algorithm (Algorithm^ for comput- 
ing certain Betti numbers (including the zero-th Betti number which is the number 
of connected components) of projections of sets defined by a constant number of 
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quadratic inequalities, without having to compute a semi-algebraic description of 
the projection. More precisely, let S C R fe + m be a compact semi-algebraic set de- 
fined by Pi > Q,...,P t > 0, with P t e R[X U ...,X k ,Y!,.. . , Y m ], deg(P) < 2, 1 < 
i < £. Let 7r : M. k+m — > R m be the projection onto the last m coordinates. In what 
follows, the number of inequalities, i, used in the definition of S will be considered 
as some fixed constant. Since, ir(S) is not necessarily describable using only qua- 
dratic inequalities, the bound in Theorem ll . II does not hold for n(S) and ir(S) can 
in principle be quite complicated. Using the best known complexity estimates for 
quantifier elimination algorithms over the reals (see we get single exponential 
(in k and to) bounds on the degrees and the number of polynomials necessary to 
obtain a semi-algebraic description of ir(S). In fact, there is no known algorithm 
for computing a semi-algebraic description of tt(S) in time polynomial in k and to. 
Nevertheless, we are able to prove that for any fixed constant q > 0, the sum of 
the first q Betti numbers of tt(S) are bounded by a polynomial in k and to. More 
precisely, we obtain the following complexity bound (see Section |3j|. 

Theorem 1.2. Let S C M. k+m be a compact semi- algebraic set defined by 

Pi > 0, . . . , Pi > 0, Pi e R[X U ...,X k ,Y u ..., Y m ], deg(Pi) < 2, 1 < % < i. 

Let 7r : R k+m R m be the projection onto the last m coordinates. For any q > 
0, < q < k, 

<? 

Y^UAS)) < {k + m)°W\ 

i=0 

We also consider the problem of computing the Betti numbers of tt(S). Previ- 
ously, there was no polynomial time algorithm for computing any non-trivial topo- 
logical property of projections of sets defined by few quadratic inequalities. We 
describe a polynomial time algorithm for computing the first few Betti numbers of 
tt(S). The algorithm (Algorithm |21 in Section[7|) computes b (n(S)), . . . ,b q (n(S)). 

The complexity of the algorithm is (A; + m) . If the coefficients of the input 
polynomials are integers of bit-size bounded by r, then the bit-size of the inte- 
gers appearing in the intermediate computations and the output are bounded by 
r(k + m) 2 ° {ql) . Note that the output of the algorithm includes b Q (Tt(S)), which 
is the number of connected components of tv(S). Alternatively, one could obtain 
bo(iv(S)), . . . , b q (n(S)) by computing a semi-algebraic description of w(S) using an 
efficient quantifier elimination algorithm (such as the one described in jSJ) and then 
using the algorithm described in [H] to compute the first few Betti numbers. How- 
ever, the complexity of this method would be worse: single exponential in k and 
m. Thus, our algorithm is able to compute efficiently non-trivial topological in- 
formation about the projection, even though it does not compute a semi-algebraic 
description of that projection (it is not even known whether such a description 
could be computed in polynomial time). 

In order to obtain Algorithm [21 we give a new construction of a certain spectral 
sequence, namely the cohomological descent spectral sequence, which plays a crucial 
role in the design of the algorithm. Even though variants of this spectral sequence 
have been known for some time |1 71 1321 ITS1 1241 [2"51 134) . to our knowledge this is the 
first time it has been used in designing efficient algorithms. The new construction 
that makes this possible is formally analogous to that of the Mayer- Vietoris spectral 
sequence, which has been used several times recently in designing algorithms for 
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computing Betti numbers of semi-algebraic sets (see |H10IE1E3!)> ancl thus this new 
construction (see Proposition 15 . 31 below) might be of independent interest. 



2. Main Ideas 

There are two main ingredients behind the results in this paper. The first is the 
use of cohomological descent, a spectral sequence first introduced by Deligne |T?ll3"2"] 
in the context of sheaf cohomology. This descent spectral sequence is used to 
compute the cohomology of the target of a continuous surjection (under certain 
hypotheses only, the limit of this spectral sequence is not, in general, the homology 
of the target). The first terms of the sequence are cohomology groups of certain 
fibered products over the surjection, and this allows to bound the Betti numbers 
of the target space in terms of the Betti numbers of those fibered products. This 
estimate was first used by Gabrielov, Vorobjov and Zell in ^5] to give estimates 
on the Betti numbers of projections of semi-algebraic sets (and more generally, 
of semi-algebraic sets defined by arbitrary quantified formulas) without resorting 
to quantifier elimination. Another use of this sequence to establish upper-bounds 
can be found in |35| which contains effective estimates for the Betti numbers of 
semi-algebraic Hausdorff limits. 

The most striking feature of this spectral sequence argument is that it enables one 
to deduce properties (for instance, bounds on the Betti numbers) of the projection 
of a set without having to explicitly describe the projection. For instance, consider 
a semi-algebraic subset of R k defined by a polynomial having a constant number 
(say m) of monomials (often referred to as a fewnomial) . It is known due to classical 
results of Khovansky (see also that the Betti numbers of such sets can be 
bounded in terms of m and k independent of the degree of the polynomial. Using 
the spectral sequence argument mentioned above, it was proved in that even 
the Betti numbers of the projection of such a set can be bounded in terms of the 
number of monomials, even though it is known (see [201) that the projection itself 
might not admit a description in terms of fewnomials. 

The construction of the descent spectral sequence given in JHj involves consid- 
eration of join spaces and their nitrations and is not directly amenable for algo- 
rithmic applications. In Section [5] we give an alternate construction of a descent 
spectral sequence. When applied to surjections between open sets this spectral 
sequence converges to the cohomology of the image. The proof of this fact is for- 
mally analogous to the proof of convergence of the spectral sequence arising from 
the generalized Mayer- Vietoris sequence. This new proof allows us to identify a 
certain double complex, whose individual terms corresponds to the chain groups 
of the fibered products of the original set. The fibered product (taken a constant 
number of times) of a set defined by few quadratic inequalities is again a set of the 
same type. 

However, since there is no known algorithm for efficiently triangulating semi- 
algebraic sets (even those defined by few quadratic inequalities) we cannot directly 
use the spectral sequence to actually compute the Betti numbers of the projections. 
In order to do that we need an additional ingredient. This second main ingredient 
is a polynomial time algorithm described in [Jj for computing a complex whose 
cohomology groups are isomorphic to those of a given semi-algebraic set defined 
by a constant number of quadratic inequalities. Using this algorithm we are able 
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to construct a certain double complex, whose associated total complex is quasi- 
isomorphic to (implying having isomorphic homology groups) a suitable truncation 
of the one obtained from the cohomological descent spectral sequence mentioned 
above. This complex is of much smaller size and can be computed in polynomial 
time and is enough for computing the first q Betti numbers of the projection in 
polynomial time for any fixed constant q. 

The rest of the paper is organized as follows. In Section |31 we recall certain 
basic facts from algebraic topology including the notions of complexes, and double 
complexes of vector spaces and spectral sequences. We do not prove any results 
since all of them are quite classical and we refer the reader to appropriate references 
[151 1261 ITT] for the proofs. In Section 01 we prove the estimate on the sum of Betti 
numbers (Theorem 11.2(1 of projections of semi-algebraic sets defined by quadratic 
inequalities. In Sectional we give our new construction of the cohomological de- 
scent spectral sequence In Section we briefly describe Algorithm ^ which is used 
to compute cohomology groups of semi-algebraic sets given by quadratic inequal- 
ities. This algorithm runs in polynomial time when the number of inequalities is 
constant. We only describe the inputs, outputs and the complexity estimates of 
the algorithms, referring the reader to [J] for more details. Finally, in Section we 
describe our algorithm (Algorithm [2J for computing the first few Betti numbers of 
projections of semi-algebraic sets defined by quadratic inequalities. 



We first recall some basic facts from algebraic topology, related to double com- 
plexes, and spectral sequences associated to double complexes as well as to contin- 
uous maps between semi- algebraic sets. We refer the reader to ^3 EH for detailed 
proofs. We also fix our notations for these objects. All the facts that we need are 
well known, and we merely give a brief overview. 

3.1. Complex of Vector Spaces. A cochain complexis a sequence C* = {C l \ i £ 
Z} of Q- vector spaces together with a sequence of homomorphisms 5 l : C % — -> C l+1 
for which 5 l o 5 1+1 = for all p. 

The cohomology groups, H l (C) arc defined by, 



where B^C) = Im^- 1 ), and Z\C % ) = Ker((5 4 ). The cohomology groups, H*(C), 
are all Q-vector spaces (finite dimensional if the vector spaces C % are themselves 
finite dimensional). We will henceforth omit reference to the field of coefficients Q 
which is fixed throughout the rest of the paper. 

Given two complexes, C* = (C 1 ,^) and D m = (Z? 4 ,<9 1 ), a homomorphism of 
complexes, (f> : C* — > D m , is a sequence of linear maps <j) 1 : C % — > D % verifying 
d'oft = ( /, l + 1 o ( 5 l for alii 

In other words, the following diagram is commutative for all i. 



3. Topological Preliminaries 



H\C m ) = Z , (C m )/B t (C), 
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A homomorphism of complexes, <j) : C" 
H*(C) -> H*(D'). The homomorphism 
homomorphisms <^>* are isomorphisms. 



-> £)*, induces homomorphisms, </>* : 
is called a quasi-isomorphism if the 



3.2. Double Complexes. A double complex is a bi-graded vector space 

ijez 

with co-boundary operators d : C 1 ^ — > C l 'i +1 and <5 : C*'- 7 — > C l+1 ^ such that 
d 2 = (5 2 = c?5 + (5d = 0. We say that C*'* is a first quadrant double complex if it 
satisfies the condition that C 1 '- 7 = when ij < 0. 

Given a double complex C**, we can construct a complex Tot*(C*'*), called 
the associated total complex of C'* and defined by Tot"(C*'*) = C % '\ with 

differential V n : Tot"(C"'*) — ► Tot" +1 (C"'*) given by D n = d + S. 




3.3. Spectral Sequences. A (cohomology) spectral sequence is a sequence of bi- 
graded complexes {Epi \ i,j,r £ Z, r > a} endowed with differentials d 1 ^ : Epi — > 
^j+rj-r+i gucn fast, (d r ) 2 = for all r. Moreover, we require the existence of 
isomorphism between the complex E r+ \ and the homology of E r with respect to 
d r : 

K+i = HdriK' 3 ) = ^r 1 r 

The spectral sequence is called a /irsi quadrant spectral sequence if the initial 
complex E a lies in the first quadrant, i.e. E 1 ^ — whenever ij < 0. In that case, 
all subsequent complexes E r also lie in the first quadrant. Since the differential 
d 1 ^ maps outside of the first quadrant for r > i, the homomorphisms of a first 
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i+j=e i+j=e+i p 



Figure 1. d r : Ep 3 -+ E l r + r - 3 '- r+1 



quadrant spectral sequence d r are eventually zero, and thus the groups Ep 3 are 
all isomorphic to a fixed group Ep 3 for r large enough, and we say the spectral 
sequence is convergent. 

Given a double complex C'*, we can associate to it two spectral sequences, 
'E^ J , "E^ J (corresponding to taking row-wise or column-wise nitrations respec- 
tively). 

If the double complex lies in the first quadrant, both of these spectral sequences 
are first quadrant spectral sequence, and both converge to H* (Tot* (C'' j), meaning 
that the limit groups verify 

(3.1) »E^=H n (Tof(C'-')), 

for each n > 0. 

The first terms of these are 'E x = H 5 {C'''),'E 2 = H d H s {C'<'), and "E 1 = 
H d {C'''),"E 2 = H s H d (C'). 

Given two (first quadrant) double complexes, C*'* and C*'* , a homomorphism 
of double complexes <f> ■ C*'* — » C''' is a collection of homomorphisms, <fi 1 ' 3 : C n ' 3 — > 
C 1 ' 3 , such that the following diagrams commute. 
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A homomorphism of double complexes, : C* : * — > C* : * induces homomorphisms 
: — > -E£ ,J between the terms of the associated spectral sequences (corre- 
sponding either to the row- wise or column- wise nitrations) . 

We will need the following useful fact (see page 66, Theorem 3.4 for a proof). 

Theorem 3.1. If (j) l s ' J is an isomorphism for some s > 1 (and all then E l r ' J 
and £y J are isomorphic for all r > s. In other words, the induced homomorphism, 
<j) : Tot*(C*'*) > Tot* ((?*'*) is a quasi- isomorphism. 

4. Proof of Theorem 11.21 

The proof of Theorem 11.21 relies on the bounds from Theorem 11.11 and on the 
following theorem that appears in |19| . 

Theorem 4.1. Let X and Y be two semi-algebraic sets and f : X — > Y a semi- 
algebraic continuous surjection such that f is closed. Then for any integer n, we 
have 

(4.1) b„(Y) < b 3 (W}(X)), 

i-\-j—n 

where WUX) denotes the (i + l)-fold fibered product of X over f: 

W}(X) = {(x , ...,Xi)G X i+1 | f(x ) = ■■■ = f^)}. 

This theorem follows from the existence of a spectral sequence Ep J converging 
to H*(Y) and such that El' 3 = H 3 (Wj(X)). Since, in any spectral sequence, the 
dimensions of the terms Ep J are decreasing when i and j are fixed and r increases, 
we obtain using the definition i|3.1|) of convergence: 

b n (Y)= dim(^)< £ dim(^), 

i+j—n i+j—n 

yielding inequality l|4.1|) . 

The spectral sequence Ep , known as cohomological descent, originated with 
the work of Deligne |171 132| . in the framework of sheaf cohomology. In ^1, the 
sequence is obtained as the spectral sequence associated to the filtration of an 
infinite dimensional topological object, the join space, constructed from /. For the 
purposes of Algorithm [3 we will give a different construction of this sequence (see 
Section EJ. 
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Proof of Theorem M.Si Since S is compact, the semi-algebraic continuous surjection 
7r : S — > tt{S) is closed: applying Theorem 14. II to ir, inequality (|4.1(l yields for each 
n with < n < q, 

(4.2) b n (ir(S)) < ]T b 3 {Wl{S)). 

i-\-j=n 

Notice that W*(S) = {(x , ...,Xi,y)\ P h (x t ,y) > 0, 1 < h < £, < t < i}. Thus, 
each W* (S) C ]R( l + 1 ) fc + m i s defined by + 1) quadratic inequalities. Applying the 
bound in Theorem II. II we get that, 

(4.3) &i(Wj(S)) < ((t + l)k + m) ^ 4+1 ». 
Using inequalities l|4.2|l and (|4.3|) we get that, 

i=0 

which proves the theorem. □ 



5. COHOMOLOGICAL DESCENT 

This section is devoted a new construction of the cohomological descent spectral 
sequence (already discussed in Section In Theorem l5.6l we obtain this sequence 
as the spectral sequence associated to a double complex associated to the fibered 
powers of X, rather than through the filtration of the join space. Convergence to 
the cohomology of the target space occurs when the map / : X — > Y is locally split 
(sec definition below). By deformation, we are able to extend the result to our case 
of interest: the projection of a compact semi-algebraic set fCorollarv l5.8ll . 

We will use this construction for Algorithm |2j 

Definition 5.1. A continuous surjection / : X — » Y is called locally split if there 
exists an open covering U of Y such that for all U S U, there exists a continuous 
section a : U — > X of /, i.e. a is a continuous map such that f(a{y)) = y for all 
yeU. 

In particular, if X is an open semi-algebraic set and / : X — > Y is a projection, 
the map / is obviously locally split. This specific case is what we will use in 
Algorithm as we will reduce the projection of compact semi- algebraic sets to 
projections of open semi-algebraic sets (see Proposition 15. 7JI in order to apply the 
spectral sequence. 

Recall that for any semi-algebraic surjection / : X — > Y, we denoted by Wj(X) 
the (p + l)-fold fibered power of X over /, 

Wj(X) = {(x , ...,x p )e XP +1 I f(x ) = ■■■ = /(%)}• 

The map / induces for each p > 0, a map from W?(X) to Y, sending (xq, . . . , x p ) 
to the common value f(xo) = ■ ■ ■ = f{x p ), and abusing notations a little we will 
denote this map by / as well. 
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5.1. Singular (co-) homology. We recall here the basic definitions related to sin- 
gular (co-)homology theory directing the reader to for details. 

For any semi-algebraic set X, let C,(X) denote the complex of singular chains 
of X with boundary map denoted by d. 

Recall that C,(X) is defined as follows: For m > 0, a singular m-simplex s is a 
continuous map, s : A. m — > X, where A m is the standard m-dimensional simplex 
defined by, 

rn 

A m = {{t , ...,t m ) | U >0 7 J2^ = !}• 

C m (X) is the vector space spanned by all singular m-simplices with boundary maps 
defined as follows. As usual we first define the face maps 

fm.i • A m ► A m _|_i, 
by fm,i{(to, . . . , t m )) = (to, . . . , 0, ti+l, . . . , t m +l). 

For a singular m-simplex s we define 

m 

(5.1) 5 S = ^(-irso/ m _ M . 

i=0 

and extend 9 to C m (X) by linearity. We will denote by C'{X) the dual complex 
and by d the corresponding co-boundary map. More precisely, given <f> 6 C m (X), 
and a singular (m + l)-simplex s of X, we have 

m+l 

(5.2) £(-l)V(*°/m,<). 

i=0 

If / : X — > y is a continuous map, then it naturally induces a homomorphism 
/* : C, (X) — > C, (Y) by defining, for each singular m-simplex s : A m — > X, 
/*(s) = so/ : A m — » Y, which is a singular m-simplex of Y. We will denote 
by /* : C'{Y) — > C"(X) the dual homomorphism. More generally, suppose that 
s = ( s 0j ■ • ■ i s p) : A m — » Wj (X) is a singular m-simplex of Wj(X). Notice that 
each component, Si,0 < i < p are themselves singular m-simplices of X and that 
/*(so) = . . . = /*(s p ) are equal as singular m-simplices of Y. We will denote their 
common image by /*(s). 

We will require the notion of small simplices subordinate to an open covering 
of a topological space (see [23 )■ Assuming that / : X — > Y is locally split, let W 
be an open covering of Y on which local continuous sections exist. We denote by 
V the open covering of X given by the inverse images of elements of IA, i.e. V = 
{/^(U) U e U}. We let (%(Y) be the subcomplex of C.(Y) spanned by those 
singular simplices of Y whose images are contained in some element of the cover 
U. Similarly, we let C^(X) be the subcomplex of C,(X) spanned by the simplices 
of X with image in V, and more generally, for any integer p, C. (W?(X)) denotes 
the subcomplex of C,(Wf(X)) spanned by simplices with image contained in V p+1 
for some V G V. The corresponding dual cochain complexes will be denoted by 
Cy(Y) and Cy(Wj(X)) respectively. We will henceforth call any singular simplex 
of C^(Y) and any singular simplex of C. (WJ(X)) admissible simplices. 

The inclusion homomorphism, t. : C^{Y) C,(Y) induces a dual homomor- 
phism, l* : C'(Y) — > Cy(Y). We also have corresponding induced homomorhisms, 
l' : C'{Wf{X)) C^(wf(X)) for each p > 0. 
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Proposition 5.2. The homomorphism l' : C'{Y) -> C^{Y) (resp. C*{Wf{X)) -» 

Cy(Wj? (A)) /or eachp > 0) is a chain homotopy equivalence. In particular, we have 
H*(C^{Y)) £ ff*(C*(y)) S H*(Y) and H* (C^(Wf (X))) S ff*(C"(Wj(A))) = 

Proof. This follows from a similar result for homology, see Proposition 2.21 in 23 . 

□ 

5.2. A long exact sequence. For each p > 0, we now define a homomorphism, 

S p : C"(Wj(X)) — > C'(W P+1 (X)) 
as follows: for each i,0 < i < p, define ir Pti : W P (X) — » W?~ (X) by, 

^*p,2 («^0 ) • • * ) ^p) — (*^0 ) * • • ; ^ii • • • i 

(ir p ,i drops the i-th coordinate). 

We will denote by (7^,;)* the induced map on C.(Wf(X)) -> C.{W^ X (X)) and 
let TT* pl : C'{W p f l {X)) -> C'(W P (X)) denote the dual map. For <j> e C*(WJ(A)), 
we define <5 P by, 

P+i 

(5.3) «j" = ]T(-i)v; +M 0. 

8=0 

Note that for any open covering V of X, the map 6 P induces by restriction a map 
C^(Wf(X)) -> C v (Wj +1 (X)) which we will still denote by 

The following proposition is analogous to the exactness of the generalized Mayer- 
Vietoris sequence (cf. Lemma 1 in [H]). 

Proposition 5.3. Let f : X — > Y be a continuous, locally split surjection, where X 
andY are semi- algebraic subsets ofW 1 andW 11 respectively. LetU denote an open 
covering of Y in which continuous sections of f can be defined on every U € U, 
and V denote the open covering of X obtained by inverse image ofli under f . The 
following sequence is exact. 

— (%oo -A c v {w Q f {x)) c- v (w}(x)) 

• • • ^ Cl{W p f {X)) C' V {W P+1 {X)) ^ • • • 

Proof. We will start by treating separately the first two positions in the sequence, 
then prove exactness for p > 1. 

(A) /* : CftPO - C' V {X) is injective. 

Let U € U and let s be a simplex whose image is contained in U. If a is 
a continuous section of / defined on U, the simplex t — <7*(s) is in C,(X), 
and verifies = s. Hence, /» : C^(X) — » C"(F) is surjective, so /* is 
injective. 

(B) /*(C£(y))=ker5°. 

Let e C™ (A). Any simplex s e C^{W}{X)) is a pair (s ,si) of 
simplices in C^(A) verifying /*(so) = f*{si). We then have 5°<j)(s) = 
</>(si) — </>(so). If = /*V> for some S C^ l (F), we have for any s, 

<*V(s) = - /V(*o) = - W*(«o)) = 0, 
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since we must have /*(so) = /*(si). Thus, we have f*{Cy(Y)) C kerJ . 

Conversely, if <fi is such that S°4> = 0, this means that for any pair (sq, si) 
of simplices in C^(Jf) verifying /*(s ) = j*(si), we have <?!>(so) = 0( s i)- 
Since we just proved in part (A) that /* : C^{X) — > C^(Y) is surjective, 
any element t E CY n {Y) is of the form t — /*(s) for some s E Cj£(Jf). 
Thus, we can define ip E C^(F) by ?/>(t) = </>(s), and the condition S°cj) — 
ensures that ip is well defined since its value does not depend on the choice 
of s in the representation t = /*(s). This yields the reverse inclusion, and 
hence exactness at p = 0. 
(C) <5 p+1 o 5 p = 0: 

From the definitions of the maps 7Tp +1 ii 7r p+2 j we have that for < z < 
p + l,0<j <p + 2, 

( 5 - 4 ) Tp+aj ° ^p+mO) = <f2,i+i ° <fi,i(0) if 3 < i- 

Let <f) E q?(W^(X)). Now from the definitions of P and <^ +1 we have 
that, 



E(-l)^ +1 (^ +1 >)); 

i=0 

p+1 p+2 



i=0 3=0 
p+1 



= E 



i=0 



E (- i ) l+% ;+2 J ° <+ M w + E (- i ) l+% ;+2,, ° ^+mW 

><7<i i<j<p+2 

»<j *>j 

Now using Equation 1)5 .4jl . the previous line becomes 

= E(- 1 ) l+J ^+2 J ° + E(- 1 ) J+% ;+2, l +i ° *;+i,M)- 

i<j i>j 

Interchanging i and j in the second summand of the previous line, we get 

= o <+i, i (^) + E(~ 1 ) j+% p+2 J +i ° ^+iM)- 

i<j i<j 

Finally, replacing j + 1 by j in the second summand above, we obtain 

= E(-i) i+ ^; +2 .,. ° K+i.M + E M) 1 ^ 1 ^ ° ^m(«); 
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and isolating the terms corresponding to j = i and j = i + 1 in the first 
sum gives 

+ E (^ 1 )' IH ^2 J «; + i,,(«); 

i<j—l i<j~l 

= 0; 

(since, again, by Equation l|5.4|l . we have 7r* +24 o7r* +1 i+1 = 7r* +2 i °tt* +m ). 
(D) Im(^) D Ker(<5P+ 1 ): 

Let G Ker((5 p+1 ). In other words, for each admissible singular m-simplex 
s = (s ,...,s p+l ):A m ^Wf +2 {X) 

p+2 

(5.5) ^(-l)V((so,...,s i ,...,s p+2 ))=0. 

i=0 

For each admissible singular m-simplex s of V let s* denote denote a fixed 
admissible singular m-simplex of X such that /*(s*) = s. Such a choice 
is possible since, as we proved in part (A), /* is surjective onto C^(Y). 
Let ip G Cy(Wj(X)) be defined as follows. For an admissible singular 
m-simplex t = (to, . . . , t p ) of Wf(X) we define 

t/j(t)=<P(f*(t),t ,...,t p ). 

Now for an admissible singular m simplex t = (to, . . . of 
with 

p+i 

W) = ^(-i)v; +1) ^(t) 

i=0 

p+1 

= ^(-l) i ^((t ,...,ti I ...,t P +i)) 
p+1 

=X)(- i )V((/»(*).*o,.. .,*<,..., <p+i)) 

i=0 



Now let s denote the admissible singular m-simplex of W^ +2 (X) defined 
by s = (f*(t), to, . . . , ti, . . . , t p+ i). Now applying Equation l|5.5|l . we get 

P+2 

5^(-l)V((*o, • Sp+a)) = 0. 

i=0 

Separating the first term from the rest we obtain, 
p+i 

<j>((t , . . . , * P+ i) = X;(-l)V((/*(*),*o, ■■■Jh-- -,t P +i)) = ^(*). 
This finally proves the exactness of the sequence. □ 
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5.3. The descent double complex. Now, let D*'*(X) denote the double complex 
defined by, D p,q (X) = C q (Wf(X)) with vertical and horizontal homomorphisms 
given by d q = (—l) p d q and S respectively, where d is the singular coboundary 
operator l|5.2|) and 8 is the map defined in l|5.3f> . Also, let D p ' q {X) = if p < or 
q < 0. 



C 3 (W°(X)) 
d 

C 2 {W%X)) 
d 

C x {W Q f {X)) 

d 

C»(Wj{X)) 

"d 




C 3 {Wj(X)) 
d 

C 2 (W){X)) 
d 

C l {WUX)) 



C°(Wj(X)) 

"d 




C 3 (W 2 (X)) 
d 

C 2 (W 2 (X)) 
d 

C 1 {W]{X)) 

d 

C°(Wj(X)) 

"d 




Lemma 5.4. The families of maps d and 5 make D°'* into a double complex. 

Proof. We need to check that d 2 = 8 2 = dS + Sd = 0. We know that d 2 = d 2 = 
since C'(W P (X)) is a cochain complex for all p, and we proved that S 2 — in 
Proposition l5.3l 

Now, suppose that <f> £ C 9 (Wf(X)) and let s = (s , . . . , s p+ i) be an admissible 
singular (q + l)-simplex of Wj! +1 (X). Then, 

d(5cf>)(s) = d fc(-l)V((*o, s p+1 ))j ; 

q+l p+1 

= (-i) p X3( _1 ) <+ ^( s ° ° ■ ■ ■ ' s * ° /w> • ■ ■ . s p+i ° /«•)• 

j=0 i=0 



We also have 



/ 9+1 

« (-1) P+1 E(- 1 )^( s ° ° ■ ' ■ ' s f+i ° /w') 

9+1 p+1 
3 =0 t=0 

Thus, it follows that dS + Sd — 0, so D*'' is indeed a double complex. 



□ 



If / : X — ► Y is locally split, and if V is the corresponding open covering of 
X defined in Section I5"T1 the double complex D*'* induces by restriction a double 
complex Dy*, where = C^(Wf(X)) when p > and q > and = 

otherwise. 
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The initial terms of the two spectral sequences associated with Dy * (cf. Sec- 
tion ar e as follows. The first terms of the spectral sequence 'El' 3 are 'E\ = 
Hs(D^'(X)),'E 2 = H^H$(Dy° {X)). By the exactness of the sequence in Propo- 
sition we have that the spectral sequence 'E^ J degenerates at the 1 E 2 term as 
shown below. 



'E 1 



CUY) 



C°(Y) 



o Jo 



o jo Jo To 













and, by Proposition 



H 3 (Y) 

H 2 {Y) 

H\Y) 

H°(Y) 



The degeneration of this sequence at 'E 2 shows that 
H* (Tot* (£>•'* (X))) s H*(Y). 
The initial term " E\ of the second spectral sequence is given by, 



"E 1 



H 3 (W°(X)) 
H 2 {W° f {X)) 
H\W S {X)) 
H°{W° f {X)) 



H 3 (W}(X)) 
H 2 (W}(X)) 
H\W}{X)) 
H°(Wj(X)) 



H 3 (W 2 (X)) 



H 2 (W](X)) 
H l {W]{X)) 
H°(Wf(X)) 
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Since this spectral sequence also converges to if* (Tot* (£)*/*) (JT)), we have the 
following proposition. 

Proposition 5.5. 

H*(Tof(D^')(X)) H*(Y). 

Proposition 15. 51 now implies, 

Theorem 5.6. For any continuous semi-algebraic surjection f : X — > Y , where 
X and Y are open semi- algebraic subsets of M n and W 71 respectively (or, more 
generally, for any locally split continuous surjection f ), the spectral sequence as- 
sociated to the double complex D*'*(X) with E\ = Hd(D°'°(X)) converges to 
H*{C'{Y))^H*{Y). In particular, 

(A) = H j [W}{X)), and 

(B) Si ff*(Tot*(D*'*(X))) - H*(Y). 

Proof. By Proposition 15 . 21 we have that the component-wise homomorphisms, i*, 
induces a homomorphism of double complexes, 

L .,. . D .,. ^ 

which in turn induces an isomorphism between the E\ terms of the correspond- 
ing spectral sequences. Hence, by Theorem 18.11 we have that, if* (Tot* (£>*'*)) = 
if* (Tot* (£>*'•)). The Theorem now follows from Proposition 1531 □ 

5.4. Truncation of the double complex. If we denote by D*'*(X) the truncated 
complex defined by, 

Dtf(X) = D^(X), if 0<i + i<g+l, 
= 0, otherwise, 

then it is clear that, 

(5.6) H l {Y) = H l (Tot' {D* q ''(X))), for < i < q. 

Now suppose that X C K fc+m is a compact semi-algebraic set defined by the 
inequalities, Pi > 0, . . . , Pi > 0. Let ir denote the projection map, ir : K fe+m — > R m . 
Let e > and let X C R fc+m be the set defined by Pi + e > 0, . . . , P e + e > 0. 

Proposition 5.7. (A) For e > sufficiently small, we have 

H*(WP(X))^H*(WP(X)), for allp>0, 

and H*(ir(X)) = H*(tt{X)). 

(B) The map, 7r|^ is a locally split semi- algebraic surjection onto its image. 

Proof. When e > is small, the sets X and X are homotopy equivalent and so are 
the sets n(X) and w(X) and the fibered products W$(X) and W%{X) for all p > 
(see @]). The first part of the proposition follows from the homotopy invariancc 
property of singular cohomology groups. The second part of the proposition is clear 
once we note that X is an open subset of R fc+m : projections of open sets always 
admit local continuous sections. □ 

We can combine Theorem 15.61 and Proposition 15.71 to construct, from the pro- 
jection of a compact basic semi-algebraic set, a double complex giving rise to a 
cohomological descent spectral sequence. 
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Corollary 5.8. Let X C M, k+m be a compact semi- algebraic set defined by Pi > 
0, . . . , Pi > and it : R k+m — > R m the projection onto the last m co-ordinates. The 
spectral sequence associated to the double complex D*' m (X) with E\ = H ( i(D*'*(X)) 
converges to H*(C'(ir(X))) = H*(ir(X)). In particular, 

(A) /:;••' =H j (Wj(X)), and 

(B) E x a JT(Tof (£>"•"(*))) = H*(tt(X)). 

Remark 5.9. Note that it is not obvious how to prove directly an exact sequence 
at the level of singular (or even simplicial) cochains for the projection of a compact 
set, as we do in Proposition ^. 31 in the locally-split setting. One difficulty is the fact 
that semi-algebraic maps are not, in general, triangulable. 

Now let X be a compact semi- algebraic set defined by a constant number of 
quadratic inequalities and / a projection map. We cannot hope to compute even 
the truncated complex Z?*'*(A) since these are defined in terms of singular chain 
complexes which are infinite-dimensional. We overcome this problem by computing 
another double complex V' q ''(X), such that there exists a homomorphism of double 
complexes, ip : V*'*(X) — > D*' m (X), which induces an isomorphism between the 
'Ei terms of the spectral sequences associated to the double complexes D''*(X) 
and T>*' m (X). This implies, by virtue of Theorem 13. II that the cohomology groups 
of the associated total complexes are isomorphic, that is, 

H* (Tot* (D''* (X))) £ H*{Tot'(V' q '{X))). 

The construction of the double complex V''*(X) is described in Section 

6. Algorithmic Preliminaries 

We now recall an algorithm described in [7], where the following theorem is 
proved. 

Theorem 6.1. There exists an algorithm, which takes as input a family of poly- 
nomials {Pi, . . • , P s } C R[A"i . . . , Xk], with deg(Pi) < 2, and a number £ < k, and 
outputs a complex T>*'' . The complex Tot* (2?"'*) is quasi-isomorphic to C t (S), the 

truncated singular chain complex of S , where 

S= p| {x g R k | P(x) < 0}. 
Per 

Moreover, given a subset T" C? , with 

S 1 = p| {x G R k | P{x) < 0}. 

PEV 

the algorithm outputs both complexes £>*'* and P/'* (corresponding to the sets 
S and S' respectively) along with the matrices defining a homomorphism Q-p^-pi, 
such that : i?*(Tot*(D*'*)) ^ H*(S) -> H*(S') ^ H* (Tot' (V''')) is the 

homomorphism induced by the inclusion i : S > S' . The complexity of the algorithm 

For completeness, we formally state the input and output of the algorithm men- 
tioned in Theorem l6.ll 

We first introduce some notations which will be used to describe the input and 
output of the algorithm. Let Q — {Qi, . . . , Q s } C ■ ■ ■ , Xk] be a family of 

polynomials with deg(Qi) < 2, 1 < i < s. For each subset J C {!,..., s}, we let 
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Sj denote the semi-algebraic set denned by {Qj > | j G J}. Notice that for each 
pair I C J C {1, . . . , s}, we have an inclusion Sj C Si. 

Algorithm 1 (Build Complex). 
Input: A family of polynomials Q = {Qi, ■ ■ ■ , Q s } C K[Xi, . . . , X&] with deg(Qi) < 

2, for 1 < i < s. 
Output: 

(A) For each subset J C {1, . . . , s}, a description of a complex F', consisting 
of a basis for each term of the complex and matrices (in this basis) for the 
differentials, and 

(B) for each pair / C J C {1, . . . , s}, a homomorphism, : F* — > F'. 
The complexes, F' and the homomorphisms 4>i_j satisfy the following. 

(A) For each Jc{l,...,s}, 

(6.1) H*(F;) = H*(Sj). 

(B) For each pair I C Jc {1, . . . , s}, the following diagram commutes. 

H*(Ff) ^ hJ \ H*(F}) 



H*(S!) H*(Sj) 

Here, (0/,j)* is the homomorphism induced by <j>i t j, the vertical homo- 
morphisms are the isomorphisms from (|6.1|l . and r* is the homomorphism 
induced by restriction. 

Complexity: The complexity of the algorithm is fc 2 ° (s) . □ 

For the purposes of this paper, we need to slightly modify Algorithm ^ in order 
to be able to handle permutations of the co-ordinates. More precisely, suppose that 
a G &k is a given permutation of the co-ordinates, and for any / C {1, . . . , s}, 
let Si.u — {(a^!), . . . ,x CT (fc)) | (x\, . . . ,Xk) G Si}. Let F° a denote the complex 
computed by the algorithm corresponding to the set 57, <r- It is easy to modify Algo- 
rithm n slightly without changing the complexity estimate, such that for any fixed 
a - , the algorithm outputs, complexes F' , F' as well as the matrices corresponding 
to the induced isomorphisms, cf)° : F' — > F' a . We assume this implicitly in the 
description of Algorithm |21 in the next section. 

7. Algorithm for projections 

Let S G M' £+m be a basic semi-algebraic set defined by 

Pi > 0, . . . ,Pt > 0, Pi S R[X lt ...,X k ,Yx,...,Y m ], 

with deg(Pi) < 2, 1 < i < £. Let it : R k+m -> R m be the projection onto the last 
to coordinates. 

The algorithm will compute a double complex, T>*'*(S), such that Tot* (£>*'* (S 1 )) 
is quasi-isomorphic to the complex Tot* (£>*'• (,!))). The double complex, V' q ''(S) is 
defined as follows. 
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We introduce k(q+2) variables, which we denote by Xij, 1 < i < k,0 < j < q+1. 
For each j, < j < q + 1, we denote by, Pij the polynomial 

P%\X\ t j, . . . , Xkj, Y\, . . . , Y m ) 

(substituting Xij, . . . , Xk,j in place of X\,...,Xk in the polynomial Pi). We 
consider each Pij to be an element of KfX^o, ■ ■ ■ , Xk, q +i, Y\, . . . , Y m ]. For each 
Pi < p < q + 1, we denote by S p C IR fe W+ 2 )+ m the semi-algebraic set defined by, 

Pi,o > 0, ■ • • , Pt,o > 0, ■ ■ • , Pi, P > 0, . . . , P Up > 0. 

Note that, for each p, < p < q + 1, and each j, < j < p we have a natural 
map, 7T p j : S p -> Sp-i given by, 

Note that in the definition above, each 5, £ M. k and tt p j exchanges the coordinates 
x j Enid x p . 

We are now in a position to define T>*>'. We follow the notations introduced in 
Section© Let Q = {Qi, . . .,Qe( q+ 2)} = {Pi,o, ■ ■ -,Pe, q +i}- For < j < q + 1, we 
let L j = {l,...,(j + l)£}c{l,...,(q + 2)£}. 

V)?{X) = Fl., 0<i+j<q + l, 
= 0, otherwise, 
The vertical homomorphisms, d, in the complex T>°'° are those induced from the 
complexes F', or zero. The horizontal homomorphisms, 5 J : F J L , — > F^ j are 
defined as follows. 

For each h, < h < i + 1, Algorithm ^ produces a homomorphism, : 
F£ s — ► F|_ , corresponding to the map ^4.1^ (see remark after Algorithm^. 

i+l 

The homomorphism <5 is then defined by, 8 = ^~^(— l) ft <fo+i,fo- We have the following 
proposition. 

Proposition 7.1. TTie complex Tot* (£>*'* (5)) is quasi-isomorphic to the complex 
Tof (£>;-• (5)). 

Proof. It follows immediately from Theorem IbM I that the columns of the complexes 
2?*'*(5) and Z?* ,, (S') are quasi-isomorphic. Moreover, it is easy to see that the quasi- 
isomorphisms induce an isomorphism between the " E\ term of their associated 
spectral sequences. Now by Theorem 13 . 1 1 this implies that Tot' (T>*'* (S)) is quasi- 
isomorphic to the complex Tot* (£>*'* (S")). □ 

ALGORITHM 2 (Computing the first q Betti Numbers). 
Input: A S C R fc+m be a basic semi-algebraic set defined by 

Pi > 0, . . . , P t > 0, 

with Pi e R[X U ...,X k ,Y u ...,Y m ], deg(Pi) < 2, 1 < i < I 
Output: b (n(S)), . . . , b q (n(S)), where tt : R fe+m -» R m be the projection onto the 

last m coordinates. 
Procedure: 

Step 1: Using Algorithm ^compute the truncated complex T>*'*(S). 

Step 2: Compute using linear algebra, the dimensions of 7^ ^ (Tot*(2?*' , )), < i < q. 

Step 3: For each i, < i < q, output, bi(ir(S)) = dim(iT (Tot* (£>••*))). 
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Complexity Analysis: The calls to Algorithm has input consisting of (q + 1)1 
polynomials in qk + m variables. Using the complexity bound of Algorithm ^ wc 
see that the complexity of Algorithm is bounded by [k + m) 2 ° lq ' □ 
Proof of Correctness: The correctness of the algorithm is a consequence of Propo- 
sition E3] and Theorem 13. II □ 



8. Conclusion and Open Problems 

For any fixed q and £, we have proved a polynomial bound on the sum of the 
first q Betti numbers of the projection of a bounded, basic closed semi-algebraic 
set defined by I quadratic inequalities. We have also described a polynomial time 
algorithm to compute the first q Betti numbers of the image of such a projection. 

Since it is not known whether quantifier elimination can be performed efficiently 
for sets defined by a fixed number of quadratic inequalities, many questions are left 
open. 

Our bounds become progressively worse as q increases, becoming exponential 
in the dimension as q approaches k. However, we do not have any examples (of 
projections of semi-algebraic sets defined by quadratic inequalities) where the higher 
Betti numbers behave exponentially in the dimension. This leaves open the problem 
of either constructing such examples, or removing the dependence on q from our 
bounds. 

Another interesting open problem is to improve the complexity of Algorithm |2 
from (k + m) 2 ° 9 to (k + m)°^ ql K Note that this would imply an algorithm with 
complexity k°^ for computing the first q Betti numbers of a semi-algebraic set 
defined by I quadratic inequalities in M. k . The best known algorithm for computing 
all the Betti numbers of such sets has complexity k 2 ° ( ' [?]• The only topological 
invariants of such sets that we currently know how to compute in time k°^> are 
testing for emptiness |2J|^ and the Euler-Poincare characteristic [0]. 
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